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Electron cyclotron maser emission mode coupling to the z-mode on a longitudinal 
density gradient in the context of solar type III bursts 
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A beam of super-thermal, hot electrons was injected into maxwellian plasma with a density gradi- 
ent along a magnetic field line. 1.5D particle-in-cell simulations were carried out which established 
that the EM emission is produced by the perpendicular component of the beam injection momen- 
tum. The beam has a positive slope in the distribution function in perpendicular momentum phase 
space, which is the characteristic feature of a cyclotron maser. The cyclotron maser in the overdense 
plasma generates emission at the electron cyclotron frequency. The frequencies of generated waves 
were too low to propagate away from the injection region, hence the wavelet transform shows a pul- 
sating wave generation and decay process. The intensity pulsation frequency is twice the relativistic 
cyclotron frequency. Eventually, a stable wave packet formed and could mode couple on the density 
gradient to reach frequencies of the order of the plasma frequency, that allowed for propagation. 
The emitted wave is likely to be a z-mode wave. The total electromagnetic energy generated is of 
the order of 0.1% of the initial beam kinetic energy. The proposed mechanism is of relevance to solar 
type III radio bursts, as well as other situations, when the injected electron beam has a non-zero 
perpendicular momentum, e.g. magnetron. 
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Electron beam injection is common in many plasma sit- 
uations. Examples are solar flares [l|, cavity magnetrons 
that generate microwaves 0, 01 or various radar applica- 
tions. Beams of relativistic, hot electrons are also likely 
to be the source of solar type III radio bursts. In the com- 
mon theory their parallel momentum is causing Lang- 
muir wave growth through the 'bump-on-tail' instability, 
and subsequent electromagnetic (EM) wave generation 
(and emission) via three main possibilities 1) non-linear 
wave- wave interactions 2) linear mode conversion on 
density gradients and 3) the antenna mechanism Q. 
Recently, particle-in-cell simulations have been carried 
out on this topic @, Q. The above mentioned mecha- 
nisms contain contribution of the parallel component of 
the beam momentum, however, in situations, where the 
beam is not completely aligned with the magnetic field, 
there is likely to be a contribution from the perpendicu- 
lar component. It is conceivable that non- field aligned, 
unstable electron beam distributions are possible in the 
solar atmosphere under some conditions. In this case, a 
bump of the distribution function will form in the per- 
pendicular direction of momentum space. If the bump is 
strong enough to achieve a positive slope, the cyclotron 
maser is triggered. The cyclotron maser is known to gen- 
erate EM emission in x-,o-, and z-mode, depending on 
the ratio of plasma frequency to electron gyrofrequency 
Q . It has been studied by PIC simulations in the context 
of the auroral kilometric radiation in Ref. [lfj • 
A beam of fast electrons is injected perpendicularly to 
the magnetic field into a 1.5D plasma, in order to turn off 
all contributions from Langmuir waves and analyse the 
purely EM signal. Time-distance plots as well as time 
evolution of emission wavelet transforms and distribu- 
tion function is investigated. A 1.5D maxwellian plasma 
is considered, allowing spatial variation in x only, while 



EM fields and particle momenta have all three compo- 
nents. In the solar atmosphere, typical values for the 
important ratio are 0.1 — 10~ 3 . We choose our pa- 
rameters as to satisfy this requirement. The background 
magnetic field is constant B = B x = 0.0003T= 3G, set- 
ting the electron gyrofrequency to uj ce = 5.28 x 10 7 Hz rad 
everywhere. We ignore the radial decrease of the mag- 
netic field in first approximation, while keeping a density 
gradient, which we believe to be crucial for the mode 
coupling mechanism here. Moreover, the constraint of 
V - B = in 1.5D does not allow for B = B(x) given that 
y— and z— components are ignorable. The background 
temperature is T — 3 x 10 5 K and isotropic. The back- 
ground plasma density at x = is no = 10 14 m -3 , giving 
u ve = 5.64 x 10 8 Hz rad. This sets ^ = 0.0935 < 1. Fur- 

ther, the electron Debye length is Xue = 3.78 x 10 _3 m, 
while the grid size is Arj e /2. The simulation setup is such 
that we study a single magnetic field line connecting Sun 
and earth. It is predicted that in the limit r 3> Rq the 
plasma density n e (r) oc r - 2 [H|, therefore, we use the 
following density profile 



n e .i{x) = n 
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with x m ax being the total system length, mimicking the 
distance of the Sun to the earth, and 
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In order to use periodic boundary conditions, the region 
Xmax/2 < x < x m ax is set to an appropriate density in- 
crease. Data analysis will, therefore, focus on the region 
< x < x max /2 only. The simulation is carried out using 
EPOCH, a fully electromagnetic, relativistic PIC code. 
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A beam of fast electrons is injected perpendicularly to 
the magnetic field at simulation time t — 0. This is be- 
cause Ref.[l2| established that the parallel momentum 
(and the associated Langmuir waves) plays no role in 
the generation of EM emission. The beam carries a to- 
tal momentum of pt = p y = rn e j^ with 7 1.155. The 
beam temperature is T b = 6 x 10 6 K. The peak beam den- 
sity is nbo = 10 n m -3 , while its spatial profile is defined 
as n b (x) = n b0 exp(-[(x - x max /2h)/{x max /AQ)f) with 

density profiles 



[c/w pe ; 



FIG. 1: beam (solid) and background (dashed) density pro- 
files at t = 0. Densities are normalized to their maxima. 

the beam density maximum at x max /2b. FigfT] shows 
the density profiles for t = 0. Note that the beam 
is injected only at t = 0. In total we simulate up to 
t = ISOWpg 1 . As a result of the above defined quantities, 
at x = 0, the plasma beta is f3 = 0.0115. The mass ratio 
is rrii/m e = 1836. The distribution functions are initially 

maxwellian. With py = p x , p± = \jp\ +p|, this leads 
to /w(p) = n b exp(-[p 2 + (p ± - p bl _) 2 ]/(2m e k B T b )) for 
the beam, where ks is Boltzmann's constant and n b is a 
normalization constant. 

For a 9 — 90° beam injection angle with respect to the 
magnetic field, in a spatial ID simulation, the beam 
is trapped at its injection point. Figj2] shows time- 
distance plots for EM field components, E x and E y , 
and changes in plasma density (left column), as well as 
spatial wavelet transforms of snapshots of E y at times 
t = 1.8, 80, 150WZ, 1 (right column), referring to the first 
pulsation maximum (top), the point where the wave 
packet passes the local plasma frequency (middle), and 
free propagation (bottom). Electric field strengths are 
given in units of Lu pe cm e /e, while distance and time are 
measured in c/uj pe and QJ~ e , 



accordingly, where ui pe is 
the plasma frequency at x — 0. The wavelet software 
was provided by C. Torrence and G. Compo [I?}- The 
wavelet transforms can be interpreted in the same way 
as described in section IV of Ref. [l2| and in more detail 
in Ref. @| 

Conventional plasma emission processes rely on the gen- 



eration of Langmuir waves in order to invoke EM emis- 
sion. Such a signal would be found in the parallel com- 
ponent of the electric field E x , as shown in Ref. [8]. It is 
evident from FigJ5J that there is no electrostatic Lang- 
muir signal produced, yet clear wave structures can be 
found in the transverse components of the EM field, see 
E y panel. The wave front of those transverse signals sug- 
gest propagation of a signal with the speed of light (units 
are chosen such that a slope of « 1 corresponds to a prop- 
agation at the speed of light), therefore, the transverse 
emission is electromagnetic. The lack of the Langmuir 
signal excludes the possibility that the generation of the 
EM emission stems from 'conventional' plasma emission 
mechanisms as mentioned in the introduction. This is 
due to the fact that in 1.5D plasma emission is not pos- 
sible [8]. Contribution from the antenna mechanism is 
also excluded by the fact that there is no density cavity 
created at any time in the simulation, as shown in the 
bottom left panel of FigJ5J Fig(3] is a representative 
snapshot of the investigated quantities of the electron 
distribution function at t = 5u;~ . The time evolution 
of the wavelet transform of the transverse electric com- 
ponent E y is shown in movie 1 [13] (including lines to 
indicate the local plasma frequency, its second harmonic 
and the z-mode cut-off at lo z — —^l 



while the respective evolution of the distribution function 
is given in movie 2 [l3[ . We can see from the distribution 
function, that the requirement for cyclotron maser emis- 
sion, g^, is fulfilled at all times. This generates waves 
at the electron cyclotron frequency, which can be seen 
in the wavelet transforms. Those frequencies are too far 
below the plasma frequency and do not allow the wave 
to propagate. It, therefore, decays again. The quick gen- 
eration and decay of waves at < t < 30^"/ appears as 
a pulsation. The frequency of this pulsation is found to 
be roughly twice the relativistic cyclotron frequency. It is 
likely that an energy exchange between beam and field, as 
shown in Ref.pjj], is happening at this stage, as the dis- 
tribution function shows subtle but complex dynamics in 
this time interval. Eventually, at rs 40kT" e , a stable wave 
packet is formed that rises in frequency, while it slowly 
moves away from the injection region. This is indicative 
of a mode coupling process on the density gradient. At 
some instance, (w, fc) of the emitted cyclotron emission 
coincides with that of the z-mode. An effect similar to the 
one discussed in e.g. Figs.l and 2 of Ref. [15;]. When the 
stable wave packet is formed, the distribution function 
shifts towards momenta higher than its initial maximum 
at « 0.5775m e c and remains rather stable. At t w 70WZ, 1 
the wave packet reaches the lower cut-off frequency of the 
z-mode, at this point the distribution function shifts to 
lower momenta again. The wave packet keeps increas- 
ing its frequency until it reaches roughly the plasma fre- 
quency at the beam injection point. 
We calculate the EM field energy by use of the Poynting 
theorem, w(t) = f[e E(x, t) 2 /2 + B(x, t) 2 /2fj, ]dx. We 
take care to exclude contributions from the background 
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FIG. 2: Left column: time-distance-plots of electric field components E x (top panel) and E y (middle panel), as well as change 
in density (bottom panel). Right: wavelet transform of E v for t — 1.8w^" e (top panel), t = 80o;^ e (middle panel), t = lSOojpe 1 
(bottom panel). Note that the background for the wavelet transform was set to white colour and does not refer to maximum 
amount of emission on the sides of the plots. Further, the black lines track the local plasma frequency and its second harmonic 



magnetic field Bq . We relate the field energy to the initial 
kinetic energy of the beam, E b k e ^ m (0) = 7.89 x 10~ 3 J. The 
overall efficiency is ss 10~ 3 E^ m (0), as shown in FigH] 
We also find the wave to be of left-handed polarization 
In the case of (oj ce /ujp e <1), the cyclotron maser is 
known to generate z-mode waves, at frequencies that are 
harmonics of the cyclotron frequency @ . The z-mode is 
left-hand polarized for cj < uj pe [l6[ . In movie 1 [l3| , it is 
shown that the wave packet is stabilized just below the 
local plasma frequency, therefore, it is no surprise to find 
the emitted wave to be of left-hand polarization. 
In Ref.[8J], it was shown that a non-gyrotropic beam in- 
jection into plasma can generate EM emission. Ref.[l2j] 
has shown that only the perpendicular momentum com- 
ponent of the injected electron beam causes t he g enera- 
tion of EM emission. The main focus of Ref.[12j was a 
parametric study, i.e. how different beam injection pitch 
angles and different gradients affect the EM emission gen- 
eration. Here the main goal is to focus on the physics of 
the EM emission generation mechanism. Thus, we car- 



ried out a 1.5D PIC simulation of a beam of relativis- 
tic, hot electrons injected with only a momentum com- 
ponent perpendicular to the background magnetic field 
into a maxwellian plasma with a density gradient. The 
choice of 1.5D turns off all EM emission generated by 
the classic plasma emission mechanism for all pitch an- 
gles, as it needs at least 2 spatial dimensions to work. 
The lack of parallel beam momentum meant that there 
was no electrostatic Langmuir wave generation, yet there 
was still EM emission due to the cyclotron maser. There- 
fore, there is no need for generation of Langmuir waves in 
order to generate EM emission via non-gyrotropic beam 
injection. The cyclotron maser generates emission at the 
electron cyclotron frequency, which - in the beam injec- 
tion region of our simulation - is far below the plasma 
frequency, hence cannot propagate and decays. Genera- 
tion and decay result in a pulsating emission generation 
until a stable wave packet is formed, that mode couples 
on the density gradient to what is likely to be a z-mode, 
until it reaches cut-off frequencies that allow for propa- 
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gation. It may well be, that the quasi-linear relaxation, 
which would result in a plateau forming of the distribu- 
tion function and eventual shutting down of the emission 
mechanism, happens on the same time scales as the in- 
verse growth rate. Calculations of growth rates are not 
part of this study. We do not see a continuous generation 
of emission, because the injected beam has no temporal 
extent (as it would have in a real solar type III burst 
situation) and beam electrons are not replenished, i.e. 
injection happens only at t = 0. Further, a propagat- 
ing beam would trigger the presented mechanism along 
its trajectory and, hence, produce the characteristic type 
III burst shape in the dynamical spectrum. Naturally, 
real life is not 1.5D, and generally there will be contribu- 
tions from classic plasma emission mechanisms and the 
antenna radiation (if density cavities are formed) as well 
as the cyclotron maser. However, in cases of strong mag- 
netic fields and straight field lines, situations can be well 
approximated by a 1.5D model. In any case, there is 
likely a perpendicular beam momentum component, that 
will contribute to the overall emission. The interplay be- 
tween EM emission generation mechanisms needs to be 



studied for various parameter spaces, but is not part of 
this work. Since the proposed mechanism has many ap- 
plications (e.g. magnetrons, radar), we believe it could 
be interesting for the physics community at large. 
The authors arc financially supported by the HEFCE- 
funded SEPNET. D.T.'s research is supported by The 
Leverhulme Trust Grant RPG-311 and STFC Grants 
ST/J001546/1 and ST/H008799/1. 

note added in proofs: After this work was complete, 
we became aware of the following: (i) when a case of 
^background = 3 x 10 8 is studied, eliminating the electron 
cyclotron maser (ECM) instability as < 0, Vp, EM 
emission identical to the one in this study still occurs; 
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FIG. 4: Total EM field energy, w(t), normalized to the initial 
beam kinetic energy, E^ m (Q). 



(ii) when a ring-shaped (in p y and p z ) beam distribution 
is considered, implying jj_ = 0, no EM emission is gen- 
erated. These findings indicate that both effects - the 
ECM instability and EM emission from transverse cur- 
rents - are present (and competing). It is probable that 
the ECM instability growth rate is so small that it can- 
not develop by the end of the simulation. Calculation 
of the growth rates commensurate to physical parame- 
ters of type III bursts and the ring distribution will be 
published elsewhere. 
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